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ABSTRACT 

Analytic models of communication networks, combining the desired 
degree of simplicity and accuracy, are noticeable by their absence 
from the literature. Modern networks are so large that the combina- 
torial possibilities impose an overwhelming burden on mathematical 
model formulation. High speed digital simulations provide one of the 
only means of solution for the systems planner. 

Military communication networks offer more hope for analytical 
solution. Much smaller than their civilian counterparts, most military 
systems can be modelled satisfactorily as aggregations of queues. 

This paper provides approximations to the probability of blocking in 
2 terminal, k-channel systems. The probability of blocking is defined 
as the probability that a vacant channel is available for connecting 
a caller to his desired party whenever the call is placed. 
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CHAPTER I 



FORMULATION OF THE PROBLEM 



1 • Introduction . 

Communication systems represent an integral part of progress in 
the social, business and technical fields. Certainly this statement 
needs no qualification when applied to the military. With the increase 
in destructive power and improvement of accuracy in modern weapons, 
large-scale dispersion of troops on the battle field has become part 
of the doctrine of warfare. Wider separation has imposed difficulty 
in controlling and directing forces at all levels of command from 
squad leader to army commander. Rapid, secure and reliable communica- 
tion is indispensable in marshalling resources to obtain the military 
and political objective. Indeed, quality and quantity of communica- 
tions systems often provides the margin of victory. 

Military communication systems encompass a wide range of sophis- 
tication; the spectrum runs from the hand-held, short range, single 
channel radio to the most modern orbiting satellite system. By far 
the most important in terms of volume of traffic and ease of control 
is the multi-channel, wide-band radio system in the Army division. 

This system closely parallels the long-haul, inter city trunking 
system operated nationally by the Bell System and independent telephone 
companies. Both of these systems have a network of radio links inter- 
connecting nodal or terminal points. 

In the Army division system, subscribers are tied by 2-wire cables 
to the switching center at the terminal. Local calls are switched and 
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connected to the called party at the local terminal. Calls to other 
terminals are transmitted over radio links to the desired terminal 
and are then connected to the 2-wire cable of the desired party. 

Civilian networks are operated with automatic switching equipment 
to route calls while military systems are at present using manual 

v 

methods in tactical operations. Another difference between the two 
systems is in the distribution of calls. The civilian system has the 
bulk of calls contained within the local area whereas the military 
network has the vast majority of calls going outside the local terminal. 
This is because subscribers in the military system are located in close 
proximity to the other subscribers in the same terminal. 

Much has been done on various aspects of communication networks. 
Perhaps the foremost source of analysis associated with developments • 
in this field is the Bell System Technical Journals. In these analyses 
efforts have been made to arrive at mathematical models admitting an 
accurate display of the significant features of communications systems. 
Practical results have been disappointing. The sheer size of a civil- 
ian system is staggering. Millions of subscribers tied into terminals 
that are connected by radio links of up to a thousand channels create 
immense combinatorial problems. Consequently, most work of a practical 
nature in the analysis of the operation of networks is done by computer 
simulation. Ease of solution provides ample compensation for the 
inelegance and inefficiency of the method. Military systems are more 
amenable to analytic solution because, as will be seen in Section 1-3, 
radio links usually have a maximum of 24 channels and rarely does a 
terminal interconnect with more than two others. Accordingly, this 
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paper will attempt to analytically model the physical operation of a 
military communications system. 

A useful parameter of a communications net is the probability of 
blocking, P^. This is defined as the probability that a call cannot be 
routed so that the desired party is reached, because all available 
channels of communication are in use. The theory of Markov processes 
and queues was used to obtain the value of this parameter. 

Approximate solutions for P^ will be obtained for a 2 terminal, 
k- channel system. In addition, an intuitive approach to the solution 
of a 3 terminal system, again interconnected with k-channel links, 
will be given. The solution for the 2 terminal model will be preceded 
by a thorough investigation of a basic system. 

2 . Organization. 

As mentioned above, this paper will attempt to develop an analytic 
model for finding blocking probabilities for small communications 
systems. An orderly progression has been attempted by initially treat- 
ing a very basic model and then, by adding real world complexities, 
arriving at a useful model for a 2 terminal system. Chapter I gives 
the background and a broad overview of the communications network with- 
in an Army division. The basic model, together with some variations 
are treated in Chapter II. More complex models, arrived at by consid- 
ering the work of operators, are introduced in Chapter III. Approx- 
imate solutions for the 2 terminal system, interconnected with a 
k-channel link, are given in Chapter IV. Chapter V provides an intui- 
tive insight into 3 terminal network problems. Areas for further 
study and extension of this paper are discussed in Chapter VI. 
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3 . The Army -Division Commualcations System . 



There is no rigid structure for the division communication system. 
Many factors must be weighed and evaluated before a concrete plan for a 
communications system can be drawn up. The type of division (armor, 
airborne, mechanized, infantry), the mission of the division, the 
geographical location, the presence or absence of other forces, the 
personnel strength; all can influence its configuration and design. 
These factors notwithstanding, there are some general remarks which can 
be made about the communications support given to the division in the 
field. The division can have from nine to fifteen battalions. Under 
present doctrine, the mixture of armor, mechanized infantry or infantry 
battalions is tailored to the mission and location of the division. 
Control is exercised by the Commanding General through three brigades 
whose strength (three to five battalions) also is regulated by local 
requirements. One brigade is normally held in reserve. 

The signal battalion, which is organic to the division, must 
furnish communication support for the division to Corps and Army HQS 
and to the main sub-ordinate headquarters. Also, within the division 
area of responsibility, signal centers must be set up to allow small 
units to connect into a division-wide communications network. Many 
modes of communication are employed. Hand-held radios suffice at the 
squad level; however as the size of the unit and the physical separa- 
tion increases, more sophisticated modes are required. These include 
line of sight microwave and tropospheric scatter. While many advances 
have been made in the past ten years, it is still true that within the 
division the main communications system is a twelve or twenty-four 
channel microwave radio system manually controlled by operators at 
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switchboards. Figure 1 shows how a typical division system might be 
configured. Usually every terminal has access to Division Main and one 
other terminal. This is a precaution against loss of a link due to 




FIGURE 1 

equipment failure or enemy action. Division Main, normally the location 
of the Commanding General and his staff is connected with as many other 
terminals (at least the brigades and division artillery) as equipment 
and personnel allow. 

At each terminal an operator is on duty to receive requests for 
service, coordinate with other operators in routing the call, and dis- 
connect the line when the call is complete. As shown in Figure 2, each 
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TERMINAL 





12 CHANNEL RADIO LINK 



SWITCHBOARD 



subscriber has an individual connection to the switchboard on the 
"subscriber" side of the board. At the switchboard, a light appears 
above the termination of a subscribers line when he turns a crank on 
his telephone set. When the operator responds to the request for 
service he determines with whom the subscriber wishes to talk. He then 
rings the operator at the terminal serving the requested party. This 
requires that at least one channel in the radio link between operators 
is not being used. When the distant operator has reached the requested 
party, the first operator tells his subscriber that the connection has 
been completed and he can commence his call. 
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CHAPTER II 



ANALYSIS OF A SIMPLE NETWORK 
DISREGARDING OPERATOR PATCH TIME 



1 . Basic Model. 

Before analyzing more complex models an understanding of the 
operation of a basic system is helpful. Accordingly the most trivial 
network will be examined first. A simple communications system is 
shown in r -gure 3. Here we have only one channel available for use 
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FIGURE 3 



between terminal A and B. The following assumptions will be made: 

a) Requests for service from the subscriber arrive at each 
switchboard with a Poisson distribution, mean rate a. 

b) The time taken to complete a conversation from the end of 
routing to hang-up, the u holding ?! time for the call, is 
distributed exponentially with mean q . 

c) Time between operator answering subscriber and the call 
commencing (the operator "patching" time) is assumed negli- 
gible in comparison to ■— ■ . 

a 

d) Having signalled the operator, subscribers wait until the 
operator answers. In other words, once a subscriber decides 
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to make a call, he is not discouraged by the time required 
to have the operator answer. 

The first two assumptions are reasonable and have been verified by a 
study on arrival rates and holding times in a civilian telephone 
exchange in New Jersey. 

This system can be treated as a basic queueing model with expo- 
nential inter-arrival time and one server. The literature generally 
denotes this as an M/M/1 queue. A steady-state solution for blocking 
probability can be obtained if the utilization rate, p , is less than 
unity where 



P = 



2a 



( 1 ) 



Steady-state implies that the system is in operation for a long period 
of time and consequently the probability of being in a particular state 
of the stochastic process is independent of the state of the process 
when first observed. This guarantees that expected system size 
(i.e., expected number waiting for service plus those being serviced) 
is always finite. The probability of blocking, p^ , in this system 

is the probability of the channel being busy and is equal to 

2a 



p b = p = 



( 2 ) 



The probability p_^ of (i - 1) callers waiting for service is 



= (1 - p) pi (3) 

It is worth noting that in this system, queue discipline is not 
necessarily first come, first served (FCFS) . Normally when queue 
discipline is not specified it is FCFS. Here, service is not neces- 
sarily rendered according to order of arrival. For example, the first 
arrival at switchboard A, with the channel busy, may occur before the 
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first arrival at switchboard B. But if the operator at B becomes 
aware of the termination of a call before his opposite number at A 
then the first subscriber in line at B will be served next although 
some have been waiting longer at A. This however is an internal 
organization problem and, while it is of interest to the individual 
subscriber, it does not alter the blocking probability. 

2. Two terminals connected by a k-channel link. 

Consider now a system consisting of two terminals with k channels 
available for communication between them. Figure 4 illustrates such 
a system. The assumptions made in the previous section will be 
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RADIO LINK 
FIGURE 4 

retained. Implicit in these assumptions is the understanding that the 
number of lines on the subscriber side of the board is so large that 
a) the future arrival rate, a , of calls is not affected by the number 
of calls presently initiated and b) the probability of a called party 
being busy can be assumed to be zero. 

The system, therefore, can be modeled as an M/M/k queue. That 
is, inter-arrival time exponentially distributed with mean — , 
service time exponentially distributed with mean — , and k servers. 
Define p^ as the probability that system size is n. Again note that 
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system size includes both subscribers waiting for service and being 
served. The probability p^ that system size is n, then can be expressed 
as 




n ^ k 

(4) 

n > k 



where 



1 - 



00 

l 

n=l 



n 



The probability of blocking, p^, is 



k-1 

p - 1 - I p (5) 

n=0 

for this system. 

3. Three terminals interconnected by 1 channel links. 



A representation of this system is shown in Figure 5. Again the 
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assumptions previously noted will be retained. One additional require- 
ment for this system is that a call can be routed over one link only. 
Here we denote the mean arrival rate at terminal i of calls going to 
terminal j as cu ^ . With operator patch time assumed zero and maximum 
route length of one link, the system may be modeled by three indepen- 
dent M/M/1 queues with 



a AB 
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01 BA 


(6) 
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a CA 


(7) 
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a BC 
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a CB 


(8) 



As before the order of arrival at the switchboard for any of these 
queues is not necessarily the order of service. 

4. Three terminal system with alternate routing available. 

Introduction of alternate routing gives a degree of flexibility 
to the previous model as well as adding realism. Alternate routing 
allows a call from A to B to be routed on links A-C and C-B, if they 
are free and link A-B is busy. Here assumptions made in Section 1 
above hold with the exception that if all links are busy incoming 
callers hang up. For ease in computation assume all arrival rates, 
a AB* a BA’ • • • 9 are the same and equal to a/ 2. Figure 6 

depicts all possible states of the system. A state is a unique 
position which completely describes the process at any instant of time. 

It is now possible to obtain a transition rate matrix A, which 
gives the rate of movement between states of the system. Solving the 
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system of equations 

P A = 0 (9) 

where P = (p^ • • • p^) the steady state probability 

that the system is in state i can be obtained. Figure 7 shows the 
transition matrix for this system. Note that it is possible to reduce 
the transition matrix from 14 x 14 to 6 x 6 by grouping similar 
states together in one "super-state". 

For example: states 2, 3, and 4 may be combined to create the 

super-state — "one call, one line busy". Using equation (8), a set 
of simultaneous equations can be written as follows: 



-a p x + 



a p 



II 
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IV 
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3 P V 
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VI 
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• (a + 3 ' P IV + T p v 
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= o 
= 0 
= 0 

( 10 ) 
= o 

= 0 

= 0 
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FIGURE 7 



TRANSITION MATRIX 
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As an example, let 



a 



1 

2 



. This yields the following results 



'II 



P III + P IV 



P V + P VI 



probability of no lines busy 



probability of 1 line busy 



probability of 2 lines busy 



probability of 3 lines busy 



112 

221 

52 

221 

(ID 

30 

221 

27 

221 



Oi 1 

For this system the probability of blocking, p^, with — — = — — is 



2 

p^ = prob {3 lines busy} + ■ prob {2 lines busy} 



27 

221 



20 

221 



57 

221 



.259 



( 12 ) 



For the same system with 2 channels per link connecting every pair of 
terminals, there are 34 possible states. Using symmetry and grouping 
similar states as before, a solution can be obtained from a set of 
11 simultaneous equations. It can be recognized that even this simple 
problem requires a significant amount of time to solve manually if 
computing equipment is not available. 

5 . Extensions .to multi-terminal networks 

It would be useful if generalizations of the results of the previous 
section were available for n-terminal networks interconnected with 
k-channel links. No such solution seems to be available. It appears 
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that systems of interconnected terminals must be treated individually 

and, except for the most trivial, the combinatorial aspects make 

[3] 

analytic solutions appear unpalatable. Benes states n the perfor- 

mance of a connecting network for a given level of offered traffic is 
determined largely by its configuration or structure . ..". This 
statement certainly has intuitive appeal and an example can exhibit 
its validity. Consider two networks, with the same number of terminals 
and total channels available for use, as shown in Figure 8. They 
differ only in network configuration. Nevertheless, certain states 




FIGURE 8 



will have different probabilities. For example consider the state 
"2 lines busy, 1 call 11 . For (i) there are eight ways to obtain this 
condition while in (ii) there are only six ways to accomplish the same 
condition even though the number of terminals and total available 
channels are the same. Obviously with different flow rates the systems 
would be incompatible. 
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CHAPTER III 



MODELS FOR NETWORKS WITH NON-ZERO PATCH TIME 



1 . Two terminal system with operator patch time . 

A more complex system to analyze is one similar to that shown in 
Figure 3 but in which operator patch time is considered significant. 
Again, we assume mean arrival rates of calls at both switchboards 
of a and holding time for calls with mean rate a. In addition operator 
patch time is taken as distributed exponentially with parameter y . 

Queue discipline is nominally first come, first serve and callers wait 
for service if the channel is occupied. 

Having only a one channel link between terminals makes available 
for this system a special analytical solution. Service time is the 
sum of two independent random variables distributed exponentially but 
with different means. Normally once an operator has completed a patch 
he is free to assist the next caller in line at the switchboard. In 
this system, however, as there is only one channel the operator is 
idle until the call occupying the channel is terminated. In practice 
the operators arrange the routing of the requested call on the same 
channel that the call itself will use. 

The system can be modeled as an M/G/l queue with arrival rate 
distributed exponentially with mean 2a. Service time, B(t), has a 
mean and variance 
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( 1 ) 

( 2 ) 
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[ 4 ] 

The embedded Markov chain technique of Kendall can be used to solve 
this type of queue and obtain answers to the most interesting questions 
(e.g. number of callers awaiting service, probability of having x 
callers awaiting service, expected waiting time for service) . Define 



X 



n 



n 
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U 



n 



n > 0 



u, , k > 0 
k = 



* 

G y (z) 
* 

F (s) 

y 



the number of calls awaiting service when 

nth call is completed 

the number of calls arriving at the 

switchboards during the routing and 
calling of the n + 1 st call. 

prob {U^ = k} 

traffic intensity 

z transform of y 

Laplace transform of y 



steady state probability that x^ = i 



We know 



[ 5 ] 



that {x } is a recurrent Markov chain if 
n 



* 

G (1) 
k 



k u, 



< 1 



k=0 



( 3 ) 



but , 



(1) = E[U ] = a m = p 



( 4 ) 



and so if p is less than unity, steady-state solutions exist, 
can be shown that 
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Gy (*> 
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(1 - p) (1 - z) F B(t) [2a (1 - z)] 



F B (t) ^ 2a n ■ z )] - z 



( 5 ) 



where 



F B(C) [2a (1 - 2)1 - 



[a + 2a (1 - z)] [y + 2a (1 - z)] 



( 6 ) 



in this example, and 
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y l 



y i 



G (z)/z=0 = 1 - p 

y 1 

I 

G (z)/z=0 

^i 

G (z)/z=0 , etc. 

y -f 



( 7 ) 



To illustrate the solution technique the following example is 
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and steady state conditions exist. 




Prob { more than 1 call at the switchboard) = 1 - y Q - 
= 1 - .4 - .23 = .37 



Expected system size is obtained from 
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Operator idle time per hour 




y l 



= .515 hours 
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Expected total time in system = — E(x^) 

= .275 hours 

Expected time waiting for an operator = .275 - m 

= .156 hours 

2 . Two terminal system with one way traffic flow. 

A much more difficult system to analyze is that shown in Figure 4 
with the following additional constraints. All traffic is assumed to 
be initiated at one terminal and exponential operator patch time (with 
mean rate u ) is considered. Now once the operator has completed a 
patch from A to B he can immediately begin to service the next customer 
(if there is one) at the switchboard providing at least one of the k 
lines is vacant. 

This system can be described as a two stage queue. Stage 1 is the 
operator patching operation and is an M/M/1 queue with infinite waiting 
room. Stage 2 is the holding time per call and is a G/M/n, n^ k, queue 
with no waiting permitted. 

The literature includes work done in the area of sequential array 
or tandem queues. This system, however, is unlike those treated 
in the literature in that the stages are not alike. They have a strong 
dependence on each other, unlike the usual tandem queues which involve 
movement from one stage to the next and where the only dependence 
of ’’successor" stages is upon their predecessor’s output. This system 
has the feature that while Stage 2 takes its arrival rate as the 
service output of Stage 1, the service rate in Stage 1 also depends 
upon the system size of Stage 2. 
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A solution technique for this system is as follows: 

Let r denote the number of channels busy (in either patching or 
conversation). Then for Stage 1, service rate for a customer is 



u 



u 

0 



r < k 
r = k 



( 8 ) 



.[7] 



Hunt L ' J shows that for a 2 stage system, both stages being M/M/1 
with unlimited waiting room, the probability of system size being n is 



n 



p (z) = E p: 

n i=0 



( 1 ) ( 2 ) 
P n-i 



where pf^ 
p.(R) 



(n +1) (1 - p) 2 p n 

= probability {i calls in stage j at steady state) 
= probability {i calls in an R-stage system) 



Stage 1 is recurrent (i.e., {p^^} , i >_ 0 is a recurrent Markov 

chain) if p, the traffic intensity or utilization rate, is strictly 

ct 

less than unity. Normally in an M/M/1 queue p = — , however, in 
this case service rate depends upon whether any lines are empty. We 
can therefore approximate u, the expected mean service rate of Stage 1 
to be 

u = u* prob {r < k) + 0* prob (r=k) (9) 

= u* prob {r < k) 

a constant. 



Stage 1 is recurrent if a <u. Then the following statement is true: 

If Stage 1 is a recurrent Markov chain the rate of departure from 

r q ] 

Stage 1 must equal the arrival rate a . J The distribution of the 
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inter-departure times is no longer exponential due to the change in 
service time distribution. However, for many problems of interest the 
output of Stage 1 can be approximated as Poisson without incurring 
great error. This is because in any real-life system the prob i // of 
lines busy = k is small. 

A series of simulations were run on an IBM 360/67 computer to 

[91 

verify this Poisson approximation. An internal package called GPSS 
was used to generate 3 runs of a 2 stage system similar to that 
described above. 1000 inputs were used in each run and reasonable 
values of a, y and a were selected. Data from these simulations are 
shown in Table 1. It can be seen that the inter departure times of 
the output of Stage 1, a ? , have the general characteristics of an 
exponential distribution (i.e., mean and standard deviation are the 
same) . In addition the mean rates of input and output of Stage 1 are 
approximately the same. By way of contrast, 2 runs were generated with 
parameter values selected so that the system was unstable. Runs #4 
and #5 of Table 1 show the difference in input and output ratio in 
Stage 1 when the service rate of Stage 2 is equal to input rate of 
Stage 1 (run #4) and when the service rate of Stage 2 is less than the 
arrival rate to Stage 1 (run #5) . 

In Stage 2 the probability of having 2 in the system is 




1 






i! 



is 



J=0 



( 10 ) 
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PARAMETERS 



RUN 


111 
a i U 0 


n 


C 

mean 


l 2 

stand. 

dev. 




i 


10 


2 


10 


6 


9.74 


9.79 




2 


3 


2 


10 


6 


2.61 


2.69 

i 




3 


5 


2 


10 


6 


4.45! 


4.56 




4 


5 


1 


30 


6 


5.94 


5.88 




5 


3 


1 


40 


6 


6.84 


i 

7.02 i 

1 


\ 



arrival rate arrival rate from stage 1 




STAGE 1 



STAGE 2 



NOTE: numbers in table represent time units between 

arrivals, departures or services 



RESULTS OF SIMULATION OF 2-STAGE QUEUE 



TABLE 1 
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and Prob {r=k} 



(ID 



k 

P 2 



k! 





where 




In Stage 1 the system utilization rate is 



P 



1 



a 



k 

u - u p 2 



k! 



k 

Z 



3=0 




k j 

= ot k! Z p 2 
j=0 j! 



( 12 ) 



k 

u k! Z p ? ^ - up3 

j=0 _£_ 

j! 

and pf 1 ^ = (1 - p x ) p 1 1 (13) 



Therefore for the 2 Stage system 




n 

Z 

i=0 




n - 



( 2 ) 



k 

z 

i=0 



P 



( 1 ) 

n-k+i 



P 



( 2 ) 

k-i 



n < k 



n > k 



(14) 
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and by using equations 8 through 14 we get 



n-i 



iSo (1 



- p x ) Pl x 



(n - i) ! ,E 0 



p 2 2 

IT 



,n < k 



( 2 ) 



n 



k-i 



i=0 (1 - p l ) P 1 



(n-k+i) 



(k - i) ! j=0 



p 2 3 

j! 



,n >_ k 



( 15 ) 



and 



(i - 



p l> i-0 



i k-i 
P 1 p 2 



(k-i)! j£ 0 



P 2 J 



(16) 



is the probability of blocking. 
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CHAPTER IV 



2 TERMINAL SYSTEMS WITH PRIORITIES AND UNLIMITED ACCESS 

1 . Models with unlimited access . 

A much more complex system to analyze is that in Figure 4 when 
arrivals from subscribers are allowed at both switchboards . As 
discussed in Section 2, Chapter III the system can be modeled as a 
tandem network consisting of a M/M/1 queue and a G/M/k queue with no 
waiting room. A subtle distinction should now be noted. The model is 
now susceptible to total collapse under the following circumstances. 
Having assumed Poisson distributed arrivals from subscribers at both 
switchboards, the probability of instantaneous arrivals at both switch- 
boards is zero. However, consider what would happen if the time 
between arrivals is small compared to the operator patch time. Operator 
B could have begun patching a subscriber before operator A could patch 
his subscriber to B. Thus the system would be suspended in a state 
of inaction with both operators trying to contact each other. Some 
real world networks prevent this situation from occuring by designating 
one channel for use by operators. Hence both operators would never be 
trying to get in touch with each other at the same time. Use of the 
"orderwire", as the designated line is sometimes called, is inefficient 
because only k-1 channels can be used for communication from A to B. 

In most systems, if the problem arises, it is disposed of by putting a 
higher priority on calls arriving at the switchboard from another 
operator than on calls from the subscriber side of the switchboard. 

If an operator is busy with one of his subscribers and notices the 
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other operator ringing he immediately stops service and attends to the 
priority call. 

Since the mean rate of calls at both switchboards is assumed to 

Ct Ct Ct 

be -y- , each switchboard in effect receives ^ high priority and — 
low priority calling rates. Operator patching time can no longer be 
considered an integral block but must be reduced to its two component 
phases: the work done by each operator. Without much error, their 

rates can be considered equal (u^ = u^ = u) . The system can then 

be depicted as shown in Figure 9. Note that exit from queues A and B 



SWITCHBOARD 

A 




SWITCHBOARD 

B 



FIGURE 9 



at 1 and 2 respectively is allowed only if the number being 
served at 3 , 4 and 5 plus the number waiting at 3 

and 4 is less than k. Queue C has no waiting (as indicated by k 
servers) . Once entry is obtained to queue C service begins immedi- 
ately. This corresponds to the fact that a call is not routed unless 
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an empty channel is available. Once the routing is complete the call 
can begin (queue C) without delay. 

2 . An approximation of the system . 

Analysis of this system is extremely complex. It combines in 
varying degree such features as tandem queues with restricted waiting 
facilities, networks of queues, and priority queues with general 
arrival distributions. All of these topics have been discussed individ- 
ually in the literature and solutions for some of the more useful 
properties obtained . 

The following approximation would be helpful in obtaining blocking 
probabilities for this system: Assume that the inter departure time 

from the queues for patching is distributed exponentially. It is 
obvious that in steady state the output of the first queue (operator) 
is no longer Poisson. Although arrivals from the subscriber side enter 
at a Poisson rate, because they have lower priority and because service 

is sometimes interrupted when k lines are full, they leave at a mean 

a 

rate of ^ but with unknown, general distribution. If, however, the 
total arrival rate a is small compared to u (effective service rate, 
see Section 2, Chapter III) the distribution of outputs entering the 
second switchboard queue should be approximately Poisson. 

To show that this approximation can be very reasonable, a queue 
having both low priority and high priority arrivals was simulated. 

Both types of arrivals were Poisson with mean rate 15 per hour. Service 
rate was exponential with mean time of 1 minute. The output distri- 
bution is shown in Figure 10. Because of the low utilization rate, 0.5, 
the inter departure time of the low priority calls has the shape of an 
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FREQUENCY OF OCCURENCE 




OUTPUT OF NON-PRIORITY 
CALLS FROM A QUEUE WITH 
SIMILAR RATE OF PRIORITY 
CALL 



MEAN 

STANDARD DEV. 



= 3.6 

= 4.1 

= 15 

= 60 




20 



TIME 



FIGURE 10 
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exponential . 

If this approximation is used the system can be modeled as in 
Figure 11. It is clear that only one operator can be busy handling 
calls from the line side of the switchboard at any one time. Assuming 
that arrivals at the switchboard from the line side are Poisson 




STAGE 1 



STAGE 2 



FIGURE 11 



distributed, entry into the waiting line for Stage 1 (call transfered 
from one operator to another) is permitted only if system size (number 
being served in Stage 2 plus number waiting in Stage 1 plus 1 or 0 
being served in Stage 1) is less than k. 

3 . System State Diagram . 

Define p_ to be the steady-state probability of having i calls 
waiting to be patched or being patched by the operator and j calls in 
progress. The state diagram for this system is shown in Figure 12 
(note that i + j cannot be greater than k) . The probability of r 
lines being busy (where 0 <_ r <_ k) is the sum of the state proba- 
bilities of the diagonal where i + j = r. 
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FIGURE 12 
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4 . System Transition Matrix* 



A precise, general analytical solution of this system does not 
appear possible. This will be more apparent later. The state diagram 
can be broken into k + 1 M super-states M , each super-state being the 
probability of r lines, 0 <_ r <_ k, being busy. While all states 
within a super-state have the same probability, aAt ,of moving to the 
next higher super-state in At, each has a different rate of movement 
to a lower state. Transition within states of a super-state is at 
equal rates u. Thus the probability of going from super-state r to 
r - 1 is dependent on which particular state of the super-state is the 
departure point. A transition matrix for rate of movement between 
super-states can be obtained as shown in Figure 13. The d^ terms are 
unknown constants, which could possibly depend upon a, a and u. 

The number of lines occupied is a "birth and death process"^^ ; 
movement is permitted only to adjacent states of the stochastic 
process at any instant of time. An expression for tt^, the steady- 
state probability of having r lines busy regardless of the state of the 
system when first observed is 



.. . a n 

r-1 r-2 



a 



0 



1 <_ r < k 



( 1 ) 



IT 



r 





or 



7T 



r 



1 



f 



r 




7T 



0 



( 2 ) 



where 



f 



r 
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SUPER-STATE TRANSITION MATRIX 



FIGURE 13 



37 



5. Determination of f . 

r 

The matrix above was solved by computer for 6 sets of parameters. 

From the results obtained for the vector tt , values for f were 

r 

obtained by noting (2) can be written 




Figure 14 shows a plot of f against r for these solutions. 

6 . Error introduced by approximations. 

To test the solutions obtained with equation 4-2 by using the 
approximations mentioned previously (i.e., assuming blocking does 
not seriously affect the distribution of inter departure times and 
disregarding the effects of priority arrivals), systems were simulated 
using the same parameters. Figures 15 through 17 compare the exact 
(simulated) probabilities of channels being used with the approximate 
solutions. In addition, a sensitivity analysis was performed on the 
approximate solutions to test the variation in blocking probability 
with changes in a. These results are shown' in Figure 18. The impli- 
cations of these results are far-reaching. It appears that the 
probability of blocking is a linear function of arrival rate, a. 

The upper boundary on a would be the minimum of patching rate or k 
times the calling rate. This information would be useful to the 
systems planner who is confronted with fluctuations in arrival rate 



over time. 
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FIGURE 18 
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CHAPTER V 



3 TERMINAL NETWORKS 



1. Simulation results . 

A simulation of a 3-terminal network was run to determine the 
effects on blocking probability of allowing routing from terminal x 
through terminal y to terminal z when all k lines between x and z are 
busy. Intuition would indicate that the mean number of channels busy 
for each pair of the 3-terminal network would be higher than that for 
a 2-terminal network with the same parameters. This is due to the 
inefficiency of routing (i.e., using 2 channels in tandem to route one 
call when the direct route is blocked) in the 3-terminal system. 
However, it could be expected that the probability of blocking for any 
call would be substantially reduced because at least 2 links must be 
at full capacity to block a call. This would happen despite higher 
individual link blocking probabilities. Table 2 lists the tt vector 
for each link of the 3 terminal system and that of a 2 terminal system 
with the same parameters. Results from the simulation showed that 
78% of calls in 3-terminal networks which were initially blocked could 
be indirectly routed to their destination immediately. 

2 . n - terminal networks with alternate routing. 

It would appear that an upper limit exists for the amount of 
alternate routing that can be used in a general, n-terminal network. 
The inefficiency of the system should increase, as a function of n, 



44 



2- 


-TERMINAL 

SYSTEM 


3-TERMINAL 
LINK 1 


3-TERMINAL 
LINK 2 


3-TERMINAL 
LINK 3 




.0019 


.0014 


.0009 


.0030 


*1 


.0049 


.0128 


.0099 


.0100 


*2 


.0319 


.0381 


.0298 


.0356 


*3 


.0618 


.0648 


.0690 


.0658 


*4 


.0917 


.1109 


.1202 


.1130 


*5 


.1565 


.1406 


.1361 


.1391 


*6 


.1525 


.1401 


.1376 


.1386 


*7 


.1326 


.1396 


.1411 


.1396 


*8 


.1276 


.1194 


.1182 


.1140 


*9 


.0737 


.0965 


.0725 


.0884 


o 
1 — 1 
t= 


.0568 


.0584 


.0621 


.0618 


*11 


.0568 


.0595 


.0511 


.0591 


*12 


.0508 


.0521 


.0506 


.0512 



TABLE 2 
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until the probability of blocking 



P 



b 




( 1 ) 



(where tt, is the probability of k channels busy in link 1) 

k l 

is greater than the probability of blocking for a 2 terminal network 
with the same parameters. Efficiency is defined as the ability of the 
network to route calls over the minimum number of links. 
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CHAPTER VI 



ANALYSIS OF RESULTS 



1. General . 

A study of Figures 15, 16, 17 indicates that the approximate 
solution for the distribution of number of lines busy in a 2-terminal 
system is satisfactory for many purposes. While the absolute differ- 
ence between the actual and approximate results is small, the differ- 
ences are significant and the approximate* solutions should be used 
for systems design or operational planning only if simulation methods 
are unavailable. The approximate results serve a purpose; they provide 
"ball-park" answers for the systems planner when formerly no infor- 
mation appears to have been available. 

The "real world" condition that after k lines are busy each 
operator can service one more subscriber was neglected in the approx- 
imate solution. This last service at the switchboard remains there 
until one of the k lines is free. Then it is patched to the other 
operator and frees the queue from which it came. This would have 
added another diagonal row of states to the state diagram, Figure 12. 

Consequently these state probabilities would be added to p. and the 

b 

approximate solution should be closer to the actual. 

The function f^ (Figure 14) appears to be very stable for the 
range of a, a and u for which the approximations are valid. Figure 
14 which shows that f^ is linear over the range of r is a very impor- 
tant discovery. From this curve, the probability of r channels busy 
can be determined by equation 4-2 if a and a are known. Extrapolation 
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for larger values of k would be possible. It is interesting to note 
the similarity between the results of the 2-terminal, k-channel 
networks with and without operator patch-time considered. Equations 
4-2 (with patch-time) and 2-3 differ only by constants. This is 
reasonable because the patch-time might be considered as just a portion 
of a larger unit of service time (patch-time and call holding-time). 
Consequently the systems differ only in mean service time. 

2. Areas for future investigation . 

Among some of the more fertile areas for further study would be 
an analytical solution for the 2 terminal, k-channel network without 
availing of the approximations used in this study. Introduction of 
priority arrivals from the subscribers into the models used here 
would also be a useful extension. It is doubtful that networks with 
more than 2 terminals would yield solutions analytically without 
benefit of over-simplifying assumptions. The nature of such assump- 
tions would also provide a fruitful investigation. 
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